Bethe Ansatz Solutions and Excitation Gap of the Attractive Bose-Hubbard Model 
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The energy gap between the ground state and the first excited state of the one-dimensional 
attractive Bose-Hubbard Hamiltonian is investigated in connection with directed polymers in random 
media. The excitation gap A is obtained by exact diagonalization of the Hamiltonian in the two- 
and three-particle sectors and also by an exact Bethe Ansatz solution in the two-particle sector. 
The dynamic exponent z is found to be 2. However, in the intermediate range of the size L where 
UL ~ 0(1), U being the attractive interaction, the effective dynamic exponent shows an anomalous 
peak reaching high values of 2.4 and 2.7 for the two- and the three-particle sectors, respectively. 
The anomalous behavior is related to a change in the sign of the first excited-state energy. In the 
two-particle sector, we use the Bethe Ansatz solution to obtain the effective dynamic exponent as a 
function of the scaling variable UL/ir. The continuum version, the attractive delta-function Bose-gas 
Hamiltonian, is integrable by the Bethe Ansatz with suitable quantum numbers, the distributions 
of which are not known in general. Quantum numbers are proposed for the first excited state and 
are confirmed numerically for an arbitrary number of particles. 



I. INTRODUCTION 

The dynamics of many simple non-equilibrium sys- 
tems are often studied through corresponding quantum 
Hamiltonians. Examples are the asymmetric XXZ chain 
Hamiltonian and the attractive Bose-Hubbard Hamilto- 
nian for the single-step growth model jjj and the directed 
polymers in random media (DPRM) [0J, respectively. 
The single-step growth model is a Kardar-Parisi-Zhang 
(KPZ) universality class growth model where the inter- 
face height h(x, t) grows in a stochastic manner under the 
condition that h(x ± l,i) — h(x,t) = ±1. The process is 
also called the asymmetric exclusion process (ASEP) in a 
different context. The evolution of the probability distri- 
bution for h(x,t) is generated by the asymmetric XXZ 
chain Hamiltonian || . The entire information about the 
dynamics is coded in the generating function (e a,l ' I,t '). 
Its time evolution, in turn, is given by the modified asym- 
metric XXZ chain Hamiltonian 



Hxxz(a) 
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with the cr's being the Pauli matrices and L the size of the 
system. (We consider only the periodic boundary con- 
dition in this work.) On the other hand, in the DPRM 
problem, the partition function Z(x, t), the partition sum 
over directed polymer configurations with fixed ends at 
(0,0) and (x,t), is the quantity of main interest. Its 
generating function (Z(x, t) n ) then evolves by the one- 
dimensional attractive Bose-Hubbard Hamiltonian [01, 



HMn) = ~\ X> 6 l+i + 6<Vi - 2) 
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Here, the 6''s(6's) are the boson creation (annihilation) 

operators, YliLi ^Pi — n is the conserved particle num- 
ber, and U{> 0) is the attractive interaction. The 
two systems are closely related, at the level of contin- 
uum stochastic differential equations, through the Cole- 
Hopf transformation, Z(x,t) = e~ h ( x '^ 0. In partic- 
ular, U in Eq. (||) is related to the particle density, 

P = Ef=i« + l)/2 of Eq. (§, by U = 4p(l - p) gj. 
Recently, the universal probability distribution function 
(PDF) at the stationary state was found for both models 
mentioned above |^|rj,[lO]| . Their common structure of the 
ground-state energies as functions of the scaling variables 
a^4Lp(l — p) and —nyLU, for the single-step growth 
model and the DPRM problem, respectively, gives the 
universal PDF at the stationary state. 

While the ground-state energy of each Hamiltonian 
gives information about the stationary state of the corre- 
sponding process, the first excited-state energy, combined 
with the ground-state energy, is related to the character- 
istic behavior of the process as it approaches the station- 
ary state. For example, in the single-step growth model, 
as t —>■ oo, the generating function takes the form 



log(e 



th{x.i) 



(3) 



(2) 



and similarly for log(Z(x,t) n ). Here, Eo(a) is the 
ground-state energy of Eq. (Q), A (a) is the inverse of the 
relaxation time, as well as the gap between the ground- 
state energy and the first excited-state energy, Ei(a), 
such that A(a) = Ei(a) — Eo(a). The size dependence 
of A (a), A (a) ~ L~ z , defines the dynamic exponent z. 

Because the asymmetric XXZ chain Hamiltonian, 
Hxxz(ce), is integrable by the Bethe Ansatz, the low- 
lying state energies, as well as the size dependence 
of the excitat ion gap, are well understood. When 
ay/ALp(\ — p) 3> 1 and the density of particles is fi- 
nite in the limit L —* oo, A(a) behaves as A(a) ~ L . 
However, when a^J \Lp(\ — p) ~ 0(1), A(a) behaves as 
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A(a) ~ L~ 3/2 |]Q. The dynamic exponent z = 3/2 
is a characteristic of the dynamic universality class of 
the KPZ-type surface growth. When the number of par- 
ticles is finite and the density of particles is very low, 
it is known that A(a) ~ L~ 2 |l2]. However, when 
a < 0, which corresponds to the ferromagnetic phase, 
most Bethe Ansatz solutions are not available although 
the Bethe Ansatz equations continue to hold. As a be- 
comes negative, the quasi-particle momenta appearing in 
the Bethe Ansatz equations become complex, so solutions 
are difficult to obtain analytically. 

The attractive Bose-Hubbard Hamiltonian is expected 
to have some resemblance to the ferromagnetic phase 
of the asymmetric XXZ chain Hamiltonian consider- 
ing the equivalence of a and —n. The equivalence is 
identified indirectly by comparing the two scaling vari- 
ables ay/ 4Lp(l — p) and — ny/ LU under the relation U = 
4p(l — p) or the two generating functions (exp(ah(x, t)} 
and (Z(x, t) n ) under the relation Z(x,t) = e~ h ^ x,t '. In 
contrast to the asymmetric XXZ chain Hamiltonian, the 
Bose-Hubbard Hamiltonian does not satisfy the Bethe 
Ansatz except in the two-particle sector JDS. Instead, 
the attractive delta-function Bose-gas Hamiltonian, 

which is the continuum version of the attractive Bose- 
Hubbard Hamiltonian, is known to be integrable by the 
Bethe Ansatz. The attractive delta- function Bose gas 
has been studied in Refs. Q and The ground-state 
energy is obtained from the Bethe Ansatz solution by us- 
ing the symmetric distribution of the purely imaginary 
quasi-particle momenta. However, the structure of the 
energy spectra is not well known for the same reason as 
in the asymmetric XXZ chain Hamiltonian with a < 0. 
The unknown energy spectra itself prevents one from un- 
derstanding the dynamics of DPRM near the stationary 
state. 

In this paper, we discuss in Section II the distribu- 
tion of the quantum numbers appearing in the Bethe 
Ansatz equation for the first excited state of the attrac- 
tive delta-function Bose-gas Hamiltonian, the knowledge 
of which is essential for solving the Bethe Ansatz equa- 
tion. In Section III, the excitation gap of the attractive 
Bose-Hubbard Hamiltonian with a small number of par- 
ticles is investigated through the exact diagonalization 
method. We show that the gap decays as A ~ L~ 2 , i.e., 
z = 2, but that the exponent becomes anomalous when 
U ~ L . The emergence of the anomalous exponent is 
explained in connection with the transition of the first 
excited state from a positive energy state to a negative 
energy state. The Bethe Ansatz solutions in the two- 
particle sector show how the behavior of the gap varies 
with the interaction. We give a summary and discussion 
in Section IV. 



II. QUANTUM NUMBER DISTRIBUTION FOR 
THE FIRST EXCITED STATE 

In this section, we study the Bethe Ansatz solutions for 
the ground state and the first excited state of the attrac- 
tive delta- function Bose-gas Hamiltonian. The eigenstate 
of H D (n) , Eq. (Q) , is of the form 

<f>(xi,x 2 , ...,x n ) 

= ^ A(P) exp(i/cpixi + ik P2 x 2 + h ik Pn x n ), (5) 

p 

where P is a permutation of 1, 2, . . . , n and x\ < x 2 < 
. . . < x n with no three x's being equal. The quasi- 
particle momenta kj's are determined by solving the 
Bethe Ansatz equations, 

k j L = 2-Kl j + J2^( 1 ^fr) (J = l,2,...,n), 

6{x) = -2 tan -1 (a:). (6) 

If the distribution of quantum numbers {/} is given, 
the set of quasi-particle momenta {k} is uniquely deter- 
mined. With such kj's, the energy eigenvalue is simply 
given by ^ = E"=i(fc]/2). 

For the ground state, the set of quantum numbers {/} 

is 

Ij = - 1 ^+j, (j = 1,2,. ..,«), (7) 

and the quasi-particle momenta are distributed symmet- 
rically on the imaginary axis in the complex-fc plane. 
Care should be taken when dealing with the first excited 
state. For the repulsive delta-function Bose-gas Hamilto- 
nian, where U is replaced by — U in Eq. (|4|), the quantum 
numbers for one of the first excited states are 

n + 1 

Ij = —+J (j = l,2,...,n-I), 



However, for the attractive case, by following the move- 
ment of the momenta as U changes sign, we find that 
the quantum numbers for the first excited state should 
be given by 

Ij =- ! ^r~+3 (j = l,2,...,n-l), 

In=-^(=h)- (9) 

That is, the two quantum numbers I\ and /„ become the 
same. Such a peculiar distribution of /j's does not ap- 
pear in other Bethe Ansatz solutions such as those for the 
XXZ chain Hamiltonian or the repulsive delta-function 
Bose-gas Hamiltonian. We remark that even though the 
two I/s are the same, all k/s are distinct; otherwise, 
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the wavefunction vanishes. Such a distribution of quan- 
tum numbers is confirmed by the consistency between 
the energies obtained by diagonalizing the Bose-Hubbard 
Hamiltonian exactly and those obtained by solving the 
Bethe Ansatz equations with the above quantum num- 
bers for very weak interactions, for which the two Hamil- 
tonians possess almost the same energy spectra. 

When there is no interaction (U = 0), all quasi-particle 
momenta, fcj's, are zero for the ground state while for 
the first excited state, all the k/s are zero except the 
last one, k n = 2n/L. In the complex-fc plane, as the 
very weak repulsive interaction is turned on, the n — 1 
momenta are shifted infinitesimally from k = with 
k± < ki < • ■ ■ < k n —i, and the nth momentum is shifted 
infinitesimally to the left from k — 2n/L. All the mo- 
menta remain on the real axis. When the interaction is 
weakly attractive, the n — 1 momenta become complex 
with Im ki < Im fc 2 < ■ • ■ < Im and Re kj ~ for 
j = 1,2,..., n— 1, and the nth momentum remains on the 
real axis, but is shifted to the left. Figure |l| shows the dis- 
tribution of the quantum numbers and the quasi-particle 
momenta in the presence of a very weak attractive in- 
teraction. The quasi-particle momenta are obtained by 
solving Eq. (||). 

Knowledge of the distribution of the quantum num- 
bers is essential for solving the Bethe Ansatz equations of 
the attractive delta-function Bose-gas Hamiltonian. For 
the original attractive Bose-Hubbard Hamiltonian, the 
Bethe Ansatz solutions are the exact solutions for the 
two-particle sector only, but are good approximate so- 
lutions in other sectors provided the density is very low 
and the interaction is very weak. This is because the 
Bethe Ansatz for the Bose-Hubbard Hamiltonian fails 
once states with sites occupied by more than three parti- 
cles are included. Thus, for the sectors with three or more 
particles, the Bethe Ansatz solutions may be regarded as 
approximate eigenstates provided states with more than 
three particles at a site do notplay an important role 
in the eigenfunctions. In Ref. |L3|, it is shown that the 
error in the Bethe Ansatz due to multiply-occupied sites 
(occupied by more than three particles) is proportional 
to U 2 , where U(> 0) in Ref. |l3| corresponds to — U in 
Eq. (||). This applies to the attractive interaction case 
also. For the repulsive Bose-Hubbard Hamiltonian, the 
Bethe Ansatz is a good approximation when the density 
is low and the interaction is strong because the strong re- 
pulsion prevents many particles from occupying the same 
site For the attractive Bose-Hubbard Hamiltonian, 
the Bethe Ansatz is good when the density is low and the 
interaction is weak because a weak attraction is better for 
preventing many particles from occupying the same site 
and because the error is proportional to U 2 . 



III. POWER-LAW DEPENDENCE AND 
ANOMALOUS EXPONENT 

We are interested in the scaling limit L — > oo with the 
scaling variable n\J\JL fixed because the common struc- 
ture of the ground-state energies of the asymmetric XX Z 
chain Hamiltonian and of the attractive Bose-Hubbard 
Hamiltonian is found in this scaling limit. In this sec- 
tion, we investigate the size dependence of the excitation 
gap by using exact diagonalization of the attractive Bose- 
Hubbard Hamiltonian in the two- and the three-particle 
sectors for L up to 30. Also, for the two-particle sector, 
we solve the Bethe Ansatz equation using the result of 
the previous section for larger L. 

Figure || shows the ground-state energies and the first 
excited-state energies versus the size of the system in the 
two- and the three-particle sectors for three values of U. 
Note that for some value of U, the sign of the first excited- 
state energy changes as the size of the system increases 
while for other values of U, no such crossover is seen in 
the range of L investigated here. 

The excitation gaps versus the size of the system are 
shown in Fig. |^ and Fig. [| on a logarithmic scale. For 
U = 0.05 and 5, the nearly straight lines indicate the 
power-law behavior of the gap, A ~ L~ z , and the slopes 
of the fitted lines indicate z ~ 2.0. However, for U = 0.5, 
the asymptotic behavior shows up only after a large 
crossover region where the effective z, z c g, is anoma- 
lously large. For the two-particle sector with U = 0.5, 
z c ff is about 2.4 in the range 14 < L < 18. For the 
three-particle sector with U — 0.5, z e ff is about 2.7 in 
the corresponding range 8 < L < 12. 

Looking into the first excited-state energy in Fig. ||, one 
can see that the anomalous value of z e g appears in the 
range of L where the first excited-state energies change 
their signs. We conjecture that the wavefunction has a 
transition as the sign of the energy changes. In order to 
confirm the connection between the anomalous exponent 
and the transition of the first excited state, we solved the 
Bethe Ansatz equation to evaluate the effective exponent 
z e g near the transition point in the two-particle sector. 

In the two-particle sector, the quantum numbers 
for the ground state and the first excited state are 
{-1/2,1/2} and {1/2,1/2} (or {-1/2,-1/2}), respec- 
tively. The quasi-particle momenta are purely imagi- 
nary, i.e., fci = —in and fc 2 = in for the ground state. 
For the first excited state, as noted already in Ref. [fTi) , 
when U < (4 / L) cos(ir / L) , the two quasi-particle mo- 
menta are real, i.e., k\ = tt/L — k and fc 2 = ir/L + k, 
while when U > (4/L) cos(tt/L), k is replaced by iq such 
that fci = 7r/L — iq and ki = n / L + iq. The examples 
of these distributions are shown in Figs. || (a) and (b) 
with L — 100 for the two values of U = 0.001 and 0.1, re- 
spectively. The transition of the first excited state occurs 
when the quasi-particle momenta are imaginary. When 
U > (4/L)cos(n/L), the ground-state energy(_E ) and 
the first excited-state energy (Ei) are, respectively, given 
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by 



and 



£ = -4sinh 2 (~), 



(10) 



i?i=4sin2 (^) cosh2 (!)- 4cos2 (^) sinh 
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where k and g are real and satisfy 
' 2 sinh k + U 



kL = log 



2 sinh k — U j 



and 



<?L = log 



£7 + 2cos(f )sinhg 
U - 2cos(f )sinh 9 



(11) 
(12) 

(13) 



respectively. We now consider the scaling limit L — > oo 
with UL finite. Let U = U* = (tt/L)su with = 
2coth(7r/2) ~ 2.181. At this value of U, the first excited- 



state energy, E\, is 0, k* 
q* = q{U*) = tt/L. Here, s . 

'2s K 



= log 



= k(U*] 
satisfies 



su 



(tt/L)s k , and 



(14) 



,2s K - su , 

which gives s K ~ 1.151. When the size of the system L 
is increased by SL with U — U* , the changes of k and q, 
5k and Sq, are, from Eqs. (|l2|) and (|l3|), 



5k = 



7TS K (4s K 2 - S(7 2 ) 



4s K 2 - 



Su 2 + (4/tt)s C ; 

2 -4) 5£ 

L2 : 



1? = 



r 6L 



SL 
L2- 



(4/ 7 r)s £/ - su 2 +4 

The perturbative expansion A(L + 5L) 
z(5L/L)), under the assumption that A(L) 
the value of z e g at U*: 

z cff = 2 ^ . 



(15) 



A(L)(1 - 
L z , gives 



(16) 



Numerical solutions of Eq. (|l4| ) give z c ff — 2.401. 

On the other hand, when U <C U* , the ground-state 
energy, E , is 0(U/L), and the first excited-state en- 
ergy, Ei, is 0(l/L 2 ), which are easily obtained from 
Eqs. ( |Io| ) and (11). Therefore, the excitation gap be- 
haves as A ~ L" 2 . Also, when U* < U < 1, 
the quasi-particle momenta of the ground state are 
±isinh _1 ([//2) — 0(e~ L ), and those of the first excited 
state are n/L ± z sinh _1 (C//2) + 0(e~ L ), which lead to 
z cff = 2. 

For arbitrary U, the effective value of z c h is evaluated 
from the relation 

_ log(A(L + l)/A(L-l)) 
Zcff " bg((£-l)/(L+l)) 

by using the solutions of Eqs. (|l^) and ([l3| ) for sufficiently 
large L. As discussed above, the exponent z c g shows an 
anomalous peak near U = U* or UL/tt = su and ap- 
proaches 2.0 as UL/tt — > or oo. Figure ^ shows a plot 
of z e ff versus the scaling variable UL/tt aX L = 10000. 



IV. SUMMARY AND DISCUSSION 



As the asymmetric XXZ chain generates the dy- 
namics of the single-step growth model, the attractive 
Bose-Hubbard Hamiltonian governs the dynamics of the 
DPRM. We studied the attractive Bose-Hubbard Hamil- 
tonian and its continuum version, the attractive delta- 
function Bose-gas Hamiltonian concentrating on the be- 
havior of the excitation gap, which is related to the char- 
acteristics of DPRM relaxing into the stationary state. 
For the attractive delta- function Bose gas Hamiltonian, 
The quantum numbers for the first excited state in the 
Bethe Ansatz equation are found for the attractive delta- 
function Bose gas Hamiltonian, and the distribution of 
the quasi-particle momenta is discussed in the presence of 
a very weak attractive interaction. Our result is the start- 
ing point for a further elucidation of the Bethe Ansatz 
solutions. We show that the excitation gap depends on 
the size of the system as a power law, A ~ L~ z , and that 
the exponent z can be calculated by using an exact diag- 
onalization of the attractive Bose-Hubbard Hamiltonian 
in the two- and the three-particle sectors and by using 
the Bethe Ansatz solution in the two-particle sector. The 
exponent z is 2.0. However, for the intermediate region 
where UL ~ O(l), the effective exponent z c s shows a 
peak. 

The equivalence of the differential equations govern- 
ing the single-step growth model and DPRM implies 
some inherent equivalence in the corresponding Hamil- 
tonians. The power-law behavior of the excitation gap, 
A ~ L~ 2 , for the attractive Bose-Hubbard Hamiltonian 
with a very weak interaction is the same as that for the 
asymmetric XXZ chain Hamiltonian with a small num- 
ber of particles, which is expected considering the rela- 
tion U — 4p(l — p). The fact that the excitation gap 
behaves anomalously for U ~ L^ 1 implies the possibility 
of an anomalous dynamic exponent z for a finite scaling 
variable ny/UL. If that is the case, one may expect the 
existence of some singularity at n ~ L -1 / 2 with finite U, 
where the dynamic exponent z takes an anomalous value, 
which is similar to the singularity at a ~ L~ x / 2 for the 
asymmetric XXZ chain Hamiltonian. 
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FIG. 2. Ground-state energies and first excited-state ener- 
gies are plotted versus the size of the system L (4 < L < 30) 
for U — 0.05, 0.5, and 5 in the two- and the three-particle 
sectors. The dotted line represents E = 0. For all values of 
U and L, the ground-state energy is negative. On the other 
hand, when U = 0.5, the excited-state energy becomes nega- 
tive near L ~ 14 in the two-particle sector and L ~ 6 in the 
three-particle sector. The signs of the excited-state energies 
for U = 0.05 and 5 do not change in the range of L shown 
here. 
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FIG. 1. For the first excited state, (a) the quantum num- 



bers Ij 's are depicted in the complex-u; plane with lj = e 



2-RiI/L 

and (b) the quasi-particle momenta kj's are shown in the 
complex-fc plane. Here, the size of the system L is 20, the 
number of particles n is 10, and the attractive interaction U 
is 0.0025. The filled circle in (a) is where the two quantum 
numbers overlap. 
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FIG. 3. Log-log plot of the excitation gaps (A) versus 
the size of the system (L) in the two-particle sector. Data for 
U = 0.05 and 5 approach straight lines with slope z — 2.0, but 
those for U = 0.5 show a strong crossover before approach- 
ing the asymptotic behavior. The solid line for U = 0.5 is 
that fitted in the range 14 < L < 18, and shows an effective 
z ~ 2.4. 
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FIG. 4. Same as in Fig. ^, but for the three-particle sec- 
tor. The fitted solid line used the data for 8 < L < 12, and 
has a slope of approximately 2.7. 
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FIG. 5. Distributions of the quasi-particle momenta, kj's, 
for the ground state (filled circles) and the first excited state 
(open circles) are shown in the complex-fc plane for n = 2. 
The size of the system L is 100 and the interaction U is (a) 
0.001 and (b) 0.1. 
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FIG. 6. Effective exponent z e ff in the two-particle sector 
versus the scaling variable UL/n at L — 10000. The interac- 
tion U varies from 0.0001 to 0.001. At UL/n = su ~ 2.181, 
z cft ~ 2.401. 
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